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ABSTRACT 



The apparently huge energy budget of the gamma ray burst GRB 990123 led to 
the final collapse of the isotropic fireball model, forcing even the most skeptical to 
consider a beamed Jet emission correlated to a supernova (SN) explosion. Similarly 
. the surprising giant flare from the soft gamma repeater SGR 1806-20 that occurred on 

qq ■ 2004 December 27, may induce the crisis of the magnetar model. If the apparently huge 

£T) [ energy associated to this flare has been radiated isotropically, the magnetar should 

\Q . have consumed at once most of (if not all) the energy stored in the magnetic field. 

On the contrary we think that a thin collimated precessing jet, blazing on-axis, may 
be the source of such apparently huge bursts with a moderate output power. Here we 

in . 

discuss the possible role of the synchrotron emission and electromagnetic showering of 
PeV electron pairs from muon bundles. A jet made of muons may play a key role in 
C"| . avoiding the opacity of the SN-GRB radiation field. We propose a similar mechanism to 

O explain the emission of SGRs. In this case we also examine the possibility of a primary 

hadronic jet that would produce ultra relativistic e^ (1-10 PeV) from n — a — e or 
neutron decay. Such electron pairs would emit a few hundreds keV radiation from 
their interaction with the galactic magnetic field, as it is observed in the intense 7-ray 
flare from SGR 1806-20. A thin precessing jet (AO ~ 10~ 9 — 10~ 10 sr) from a pulsar 
may naturally explain the negligible variation of the spin frequency v = 1 / P after the 
flare (Au < 10~ 5 Hz). A correlation between SGR 1806-20, SGR 1900 +14, Cygnus 
and the AGASA excess of EeV cosmic rays (CRs) has been found; we suggest that a 
jJJ ■ robust EeV signal may be detected by AUGER or Milagro in the near future. 



1 INTRODUCTION duced in shocks of highly relativistic outflows with different 

Lorentz factors via synchrotron radiation of electron pair. 

1.1 GRBs fireball versus thin precessing jets These so-called "internal shocks" are thought to produce the 

highest energy gamma emission. The consequent interaction 

The association between gamma ray bursts (GRBs) and SNe c , nl . , , nl ,- / . 

b J v ' of the stellar ejecta with the circum-stellar medium (exter- 

has become more and more convincing in the last few years n£j ghocks) wQuld fae responsible of the mu iti- W avelength 

after the discovery of two main events. In 1998, the spec- e, , , , ,. . , , c 

■' ' v afterglow and its complex time structure discovered tor the 

troscopical and photometrical analysis of the GRB 980425 firgt ^ by BeppQ gax (Cogta ^ ^ lgQ7) 
afterglow revealed for the first time the presence of a SN 

related to the GRB event (Galama et al. 1998), and more However isotropic explosions would require extremely 

recently SN2003dh has been detected in coincidence with high input energies, the larger and harder the more distant 

the GRB 030329 (Hjorth et al. 2003; Stanek et al. 2003). the GRBs (the so called "Amati " law, whose beahaviour 

At least other two GRB events seem to be connected to is opposite to the Hubble law); for example, to account for 

SNe explosions (GRB 011121/SN2001ke, Garnavich et al. the GRB 990123 emission, about two solar masses should 

2003; GRB031203/SN20031w; Malesani et al. 2004). At the have been entirely converted in gamma-rays (which should 

time of its discovery, the GRB980425/SN1998bw association be doubled if MeV neutrino pairs were taken into account; 

faced large resistance among the scientific community, whose Fargion 1999). Very massive black holes (10 - 100 M©) are 

attention was focused on the fireball scenario where this needed to fuel such fireball explosions. But the inner GRB 

SN/GRB connection was unacceptable. Indeed the spher- time structure, as short as a tenth of milliseconds, would re- 

ical fireball has been in the last decade the most popular quire a Schwarzschild radius smaller than a few tens of km, 

model to interpret GRBs, with the basic assumption that corresponding to objects of only a few solar masses. Such a 

the energy is radiated isotropically. Consequently the out- discrepancy, added to the SN/GRB connection, made clear 

put power of GRBs would be billion times higher than that the inconsistence of the fireball scenario. Later on new fam- 

of a SN. According to the fireball model, GRBs are pro- ilies of fireball models with a wide beamed emission have 
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been introduced to overcome this puzzle and reduce the en- 
ergy budget (such as the Hyper-Nova, Supra-Nova, or the 
Collapsar models). In these "compromise" scenarios the out- 
flow was thought to occur in a "fountain" cone with opening 
angles as large as ~ 5° - 10° . However we think that even this 
"mild" beaming is unable to solve the whole GRB puzzle. 

In fact the beamed emission proposed in such models 
reduces the required power output by three orders of magni- 
tude, down to 10 50 erg s _1 , nonetheless this is yet 10 6 - 10 7 
times more intense than the observed SN power. There is no 
reason to expect such a huge unbalance between the GRB 
jet and the SN luminosities, Egrb — 10 6 x Esn- Moreover 
a few degrees jet does not solve the puzzle of the low 7 ray 
power detected in the GRB980425/SN1998bw event, which 
finds a natural understanding assuming a narrow jet blazing 
the observer off-axis (Wang & Wheeler 1998; Fargion 1999). 
Moreover a single one-shot explosion beamed in a wide jet 
cannot account for the presence of X-Ray precursors before 
the onset of a few GRB events (~ 6%). 

In this context, we have been proposing for the last 
decade, an alternative scenario to the fireball where GRBs 
(as well as SGRs) are originated by highly collimated pro- 
cessing and blazing gamma jets with an aperture angle 9 ~ 
0.02° - 0.006° (Fargion & Salis 1995a, 1995b, 1995c; Fargion 
1999). 

The much narrower solid angle of the jet, AQ/O, ~ 10 -8 
- 10~ 9 , can reconcile (assuming the energy equipartition be- 
tween the powers of the SN and the jet) the puzzle of the 
"low" observed isotropic SN luminosity respect to the ap- 
parent "huge" GRBs power (an observed 10~ 8 - 10~ 9 ratio). 
The angular dynamics (a spinning and multi-precessing jet) 
may explain the wide range of properties observed in differ- 
ent GRB events (multiple bursts, re-brightening and bumps, 
X-ray precursors, dark GRBs, X-ray flashes, low luminous 
off-axis GRBs; Fargion 1999; Fargion 2003). 

Within this framework, the Amati-Ghirlanda E pea k — 
Ei SO correlation (Amati et al. 2002, Ghirlanda et al. 2005) 
is just the effect of a very biased selection. At increasing 
distances (high z) and volumes, a larger sample of highly 
collimated GRBs becomes observable; the jets more colli- 
mated towards the observer appear as very hard and pow- 
erful GRBs, while those blazing off-axis would be mostly 
obscured and below the detectors' sensitivity. At extreme 
redshifts the number of GRBs should drop off (due to a fi- 
nite age of their progenitor and the absence of SN events) or 
they might correspond to the mysterious short GRB events. 

However assuming a thin one-shot jetted model leads 
to a contradiction with the observations and it does not 
provide the ultimate solution to the issue. In fact, such a 
strong collimation would imply that the rate of one-shot 
beamed GRBs should be F^ ed ~ 10 8 - 10 9 rg£ B , where 
Tors ~ 10~ 7 - 10~ 8 yr" 1 galaxy" 1 . Given that GRB are 
expected to be related to SN explosions, this value would be 
at least 10 2 — 10 4 times higher than all the observed SN rate, 
Fgjv ~ 10~ 2 yr _1 galaxy _1 . This apparent inconsistence may 
be solved if one assumes that the jet structure is precessing 
and persistent, with a decay lifetime much longer than the 
typical observed GRB timescale {tcRB ~ 10 s, tj e t > 10 4 s, 
Ej et — Ep ea k(jj— t )~ a , where a > 1, in order to converge 
to the value of the total output energy). Therefore we argue 
that the jet continues to exist after the SN event, although 



it is fading with time. This reduces the number of GRBs 
needed to account for the observed GRB rate. The activity 
of the jet might be later rejuvenated by an accretion disc, 
possibly fed by a nearby compact companion. The torque 
perturbation created by the companion induces the multi- 
precession of the GRB jet. 

This might also explain the abundance (nearly two- 
third) of the so-called "orphan" or "dark" GRBs, where the 
brightest event associated to the SN and its consequent af- 
terglow may be already faded just a few weeks or months 
before. The delayed jet activity may blaze again and it may 
be source of rare optical and radio re-brightening of the 
GRB, while much older and weaker sources, may be observed 
thousands of years later as SGRs, or Anomalous X-Ray pul- 
sars 1 , provided that they are nearby and on-axis. This occurs 
mostly in our galaxy and in the Large Magellanic Cloud. 

Let us remind that a few GRBs and SGRs did show 
comparable time structures and spectra (Fargion 1999), even 
if a possible connection between GRBs and SGRs has been 
often underestimated or even rejected (Woods et al. 1999). 
Consequently the "necessary" evolution of the fireball sce- 
nario into a jetted model may be reflected on SGRs, whose 
blazing nature challenges the interpretation given by the 
magnetar model (Duncan & Thompson 1992; Thompson & 
Duncan 1995). In fact, similarly to the GRB 990123 fireball 
crisis, the very exceptional giant flare from SGR 1806-20 oc- 
curred on 2004 December 27 is suggesting that an isotropic 
magnetar has to be questioned as the explanation of these 
events for several reasons that we discuss in the following 
section. 



1.2 SGRs and the puzzles of the magnetar model 

Soft Gamma Repeaters are nearby X-ray sources localised 
within the Milky Way (excluding SGR 0525-66 in the Large 
Magellanic Cloud), emitting sporadic short bursts of soft 
gamma photons up to hundreds of keV. The giant flare from 
SGR 1806-20 on 2004 December 27 has been characterised 
by an energy release which largely exceeded all previous 
recorded events. It should be noticed that the giant flare 
energy exceeded the integral sum of all GRBs and SGRs 
ever recorded. This unique event, averaged over 100 years, 
would make the SGRs flux much more intense than that of 
GRBs (< $sga >^ 100 eV cm" 2 s" 1 sr" 1 ). The flare, if it 
was blazing far from the sun in a dark night (contrary to the 
real case), would have shone (while showering and ionising 
the atmosphere) like a brief visible aurora in the sky. 

At an assumed distance of 15 kpc the isotropic energy 
associated with the flare would be E = 3.5 — 5 x 10 46 df 5 
ergs (Hurley et al. 2005; Schwartz et al. 2005) with a peak 
luminosity in the first 0.125 s equal to L = 1.8 — 2.7 x 10 47 
erg s _1 . Only two other giant outbursts of energy above 10 44 
ergs have been observed in about 30 years of previous SGR 
activity, but never at such a huge peak power (E ~ 10 47 erg 
s _1 ). While the giant flare energy is a tiny fraction (10~ 4 ) 
of an electromagnetic galactic Supernova, its peak power 
exceeds the optical SN output by nearly the inverse ratio(~ 

1 There might be a natural extension of the jet model to nearby 
low energy sources such as the X-Ray Flashes (XRF), where a jet 
blazing off-axis could explain their properties. 
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10 4 ). Note that the energy budget could be lowered if SGR 
1806-20 were closer, as it has been recently proposed that 
its distance may be within 6.4 - 9.8 kpc (Cameron et al. 
2005). Note also that a closer distance implies an extreme 
relativistic propagation of the radio ejecta (> 0.6c). 

The most popular model to interpret the properties of 
SGRs is the magnetar, a neutron star with a very high 
magnetic field (B ~ 10 14 - 10 15 G). This model implies an 
isotropic thermal energy release like a mini-fireball, and it 
was born in the 90's when fireball models were ruling the 
understanding of GRBs. Surprisingly, a flare of energy 3.5 
- 5 x 10 46 d 15 erg radiated isotropically would need to con- 
vert at once almost all of the magnetar energy stored in the 
magnetic field. This would imply a remarkable decrease of 
the magnetic field energy density which should appear as a 
prompt decrease in the pulse period times the period deriva- 
tive, PP. As a consequence, there should be a correlation 
between the burst activity and the variation of the spin pe- 
riod, that so far has not been observed (Woods et al. 2002). 
Neither the recent giant flare showed evidence for a variabil- 
ity in its period (Mereghetti et al. 2005, Woods et al. 2005 
ATEL n. 407). 

Moreover the high energetics involved in this intense 
flare has raised the possibility that nearby extragalactic 
SGRs may be connected to short GRBs. In this case if 
SGRs are isotropic phenomena, as predicted by the mag- 
netar model, one should observe an anisotropy in the dis- 
tribution of short GRBs compared to long ones, with the 
first more concentrated towards local cluster of galaxies, es- 
pecially towards Virgo the closest rich cluster in the local 
universe. Virgo contains about 128 spiral galaxies (Binggeli, 
Tamman & Sandage 1987), assuming that the periodicity of 
a giant flare in a galaxy may be 1 in 100 years, one should 
expect to see 128/100 ~ 1.3 events/yr clustered in the area 
of Virgo. If short GRBs and SGRs are the same events, the 
BATSE instrument in its 10 years life time, should have de- 
tected about 13 short GRBs with a spatial distribution cor- 
related to the Virgo cluster. Such a remarkable anisotropy, 
expected in the magnetar scenario, does not appear in the 
BATSE GRB catalogues neither from the Virgo cluster, nor 
from galaxies with extremely high star formation rate such 
as Arp 299 and NGC 3256 (Popov & Stern 2005). Moreover 
no similar giant flares were found in the BATSE data in 
nearby galaxies (within 2 and 4 Mpc) with ongoing massive 
star formation such as M 82, M 83, NGC 253 and NGC 4945 
(Popov & Stern 2005). According to the authors, the galac- 
tic rate of giant flares with energy around 10 46 erg should 
be less than 10~ 3 yr _1 . 

On the other hand, a very thin, precessing 7 Jet would 
reduce the high energy budget required in an isotropic en- 
ergy release, and it may produce at any time an exceptional 
bright (on-axis) event. The rarety of such an event would 
imply a negligible activity even from Virgo, still consistent 
with the absence of an anisotropy in the BATSE short bursts 
catalogue; the giant flare visibility may be below the thresh- 
old from clusters at larger distances (z > 0.05) (Popov & 
Stern 2005). 

1.3 Structure of the paper 

We suggest that precessing and persistent jets would better 
explain the observed properties of both GRBs and SGRs, 



even though several issues are yet to be addressed within 
this model. In this paper we aim to improve our scenario 
introduced and discussed in previous works (Fargion 1998; 
Fargion 1999; Fargion 2003), where we considered jets made 
of relativistic GeV electrons producing the gamma radiation 
by ICS, and we try to solve some of the problems related to 
this model. How can a jet of GeVs electron pair penetrate 
and propagate through the SN envelope? How can the col- 
limated pairs survive the dense photon background at the 
baryon column depth at the SN surface? Is it possible that a 
jet is already present before the SN explosion, but it reaches 
its maximum power only after the stellar collapse? 

To overcome the opacity problem due to the photon SN 
background, here we present a scenario for GRBs where the 
jet is made by muons of energy > 1 - 10 PeV. We assume 
(but we will not thoroughly discuss this hypothesis in this 
paper) that in GRBs PeV muons are themselves secondaries 
of ultra high energy (UHE) neutrinos possibly originated in 
a deeper and inner jet made of EeV - ZeV nucleons accel- 
erated in the GRB/SN explosion. We show that the PeV 
muons could more easily propagate through the dense SN 
radiation background and they could also better penetrate 
the outer stellar layers compared to a jet made of electrons 
(or baryons). When muons decay into ~ PeV electrons, wc 
compare the energy loss mechanisms that can produce the 
observed gamma emission: synchrotron interactions due to 
the stellar or galactic magnetic field, versus ICS of e ± onto 
the stellar radiation background. 

We discuss a similar mechanism to explain the gamma 
signal from SGRs, always assuming a primary jet made of 
muons. We find that synchrotron losses (due to the galac- 
tic magnetic field B ~ 2.5 x 10~ 6 G) of electron pairs at 
E e ~ PeVs, from the decay of PeV muon bundles, may be 
the best mechanism to originate the 7 radiation from SGRs. 
We consider also that PeV muons may be secondaries of 
EeV nucleons originated via photopion interactions, since we 
find a possible link between the anisotropy in the EeV CR 
distribution near the galactic center observed by AGASA 
(Hayashida et al. 1999) and SUGAR (Bellido et al. 2001), 
and the position of SGR 1806-20. If this correlation is not 
fortuitous this maybe the first evidence of a connection be- 
tween EeV CRs and the SGR activity. The possibility of 
detecting CRs (nucleons and TeV photons) and neutrinos 
from this giant flare in AUGER, MILAGRO and AMANDA 
is also discussed. 

The paper is organized as follows: in §2 we briefly sum- 
marise the main characteristics of the precessing jet model, 
in §3 we discuss the problem of the opacity of electrons prop- 
agating through a dense supernova gas photon, and we in- 
troduce a primary muon jet model, in §4 we discuss the 
possible mechanisms that can lead to the gamma emission 
from secondary electron pair Jets, in §5 we apply the same 
model to SGRs, and we discuss the correlation between EeV 
CRs and SGR 1806-20, and finally in §6 we present our final 
discussion and conclusions. 



2 THE PRECESSING JET SIGNATURES 

Precessing jets have been already introduced to explain the 
properties of other astrophysical objects. Well known micro- 
quasars, such as 1915-16 and SS 433, are the best candidate 
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cretion disc or a companion in the case of SGRs. In SGRs 
the close encounter between a compact source with a fading 
jet and the companion star, may strip matter from the lat- 
ter, that, falling into the accretion disc, leads to an increase 
or a revival of the jet power. The presence of a compact star 
(white dwarf or a neutron star) would be the best candidate 
to bend the jet before and after a SN in GRBs. 

Therefore we propose a model of a jet that is not static, 
but we assume that it is spinning as the dense central ob- 
ject that hosts it (pre-SN star or NS), and precessing, even 
at random, due to the presence of a binary companion (or 
an asymmetric accretion disc ). A nutation mode is super- 
posed to the precession circle. The spinning motion provides 
a short duration of the GRB (or SGR) event, while precess- 
ing modes move the jet in and out of the line of sight leading 
to a re-brightening or a quasi periodic activity. 

The temporal evolution of the angle between the jet 
direction and the rotational axis of the object, 9\(t), can be 
expressed as 



j0i««5« J0!«4*K6 J02««6S« JC2««MI ( 




Figure 1. Up: The Egg Nebula, whose shape might be explained 
as the conical section of a twin precessing jet interacting with the 
surrounding cloud of ejected gas. Down: The observed structure 
of the outflows from the microquasar SS433. A kinematic model 
of the time evolution of two oppositely directed precessing jets is 
overlaid on the radio contours (from Blundell & Bowler 2005). 



to show the multi-precessing jet evolution in nearby galac- 
tic objects (see Fig. 0. Other well known sources as CH 
Cyg, Cyg X-3 and GRO J1655-40 show precessing jets whose 
properties remind, at a smaller scale, those of larger active 
galactic nuclei such as M87 or Cygnus A. Bipolar planetary 
nebulae as the Sa2-237 (Schwarz, Corradi & Montez 2002) 
hidden in molecular clouds or the Egg Nebulae (see Fig. 0, 
also show features that give evidence of a twin double cone 
possibly made by a precessing jet. At a lower luminosity, 
jets are used to explain the properties of two accreting white 
dwarfs (Zamanov et al. 2005) exhibiting collimated outflows 
with a morphology well fit by a precessing jet model. Out- 
flows collimated in jets are also observed in Herbig Haro 
protostellar objects (HH 40, HH 34), and in the large pro- 
tostar IRAS 16547-4247. 

According to Mirabel (2004) there is an universal mech- 
anism for the production of relativistic jets in accreting black 
holes that would unify AGN, GRBs, microquasars and Ultra 
Luminous X-ray (ULX) sources. 

In our model we assume that GRBs and SGRs originate 
from the same process and they represent the early and the 
late stages of the evolution of a precessing jet. The jet is 
possibly fuelled by the SN event for a GRB and by an ac- 



0i (t) = y/oi + e 2 y (i) 

with 

Ox (t) = Sm(Wbt + 4>b) + Opsr ■ sin(Wp»rt + <j>psr) + S ■ 

sin(w s i + (f> s ) + On ■ sin(cjjvt + </>n) + X (0) 

y (t) = 9 a ■ sintJot + cos(iv b t + (f> b ) + 6 psr ■ cos(u) p3r t + 

<t> V sr) + 3 ■ COs(ui s t + 4> s ) + 9 N ■ COs(ui N t + <j>Jf)) + 0y(0) 

7 is the Lorentz factor of the relativistic particles (half 
PeV electrons). psr , 9 b , 9n are respectively the maximal 
opening angles due to the spinning of the star, the pertur- 
bation due to the companion, and the nutation motion of 
the multi-precessing jet axis. The arbitrary phases 4>b, cipsr, 
4>n, for the binary, spinning pulsar and nutation, are able 
to fit the complicated features of the GRB light curves in 
most of the observed events. The additional parameters have 
been introduced to model the light curve of the SGR giant 
flare. 9 a is the bending angle of the entire cone precessing 
with a frequency ui a , corresponding to a period of nearly a 
dozen years (related to "intermittent" SGRs activity) 2 ; 9 S is 
a possible inner pulsation of the pulsar at frequency u s = 25 
rad/s; 9n and uin refers to the nutation motion of the jet, 
and 6 X (0), 0j/(O) are fine tuned phases able to reasonably 
reproduce the observed light curve of the giant flare. 

In Fig |21 we show the results of our models for the pa- 
rameters displayed in Tabled 

We have shown already (Fargion 1999, Fargion 2003) 
that by varying the parameters of this simple model one 
can fit the profiles of the GRB events. Moreover we proved 
(Fargion & Salis 1995a, 1995b, 1995c, Fargion & Salis 1998, 
Fargion 1999) that the GRB 7 observed spectrum may be 
well fit by the Inverse Compton Scattering of relativistic 
electrons with a power-law or a monochromatic spectrum. 



2 The real binary angular velocity might be very slow: in this 
case LU a might play the role of the the binary angular velocity 
while the other ones {wi,, u>n) are shorter precessing and nutating 
frequencies due to the asymmetry of the spinning jet. 
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Figure 2. A possible 3D structure of the precessing jet obtained 
with the parameters displayed in Table 1; at its centre the "ex- 
plosive" SN-like event that here we have represented with the 
Eta-Carinae lobes, or a steady binary system system where an 
accretion disc around a compact object powers a colimated pre- 
cessing jet. 



7 = 10 u 8 a 

6b — 1 ®pi 

ti?i,=4.9 TO - 4 rad/s u) p 
(j> b = 2tt - 0.44 



0.2 

=1.5 -10 7 /7 
=0.83 rad/s 

psr— 7T + 7T"/4 



1.6 ■ 10" 8 rad/s 



S =1.5 -10 6 /7 



9 N =5 ■ 10 7 / 7 
wjv=1.38 10 -2 rad/s 
0jv=3.5 7r/2 + 7r/3 



: 25 rad/s 



Table 1. The parameters adopted for the jet model represented 
in Fig.H 

3 THE NATURE OF THE GRB JETS 

3.1 The opacity of the SN radiation to an 
electron pair jet 

One important issue is whether the jet is born during the 
SN event or if it is already present. 

If the jet is created during or immediately before the 
SN explosion as soon as the GRB emission is detected, it 
is unlikely that a jet powered for the typical duration of a 
long burst ~ 20 s could penetrate a giant star. It would 
rather dissipate its energy in the stellar envelope. Moreover 
it would take about 1000 s for the jet to penetrate a red 
supergiant and break the stellar surface (Rosswog 2003). 

Recent theoretical models assume that the progenitor 
of a GRB is a massive star which has already lost its H 
envelope before core collapse, given the correlation between 




Figure 3. The jet's 2D projection along a plane orthogonal to 
the vertical axis. Such a configuration of the jet has been obtained 
with the parameters displayed in Table 1. 



GRBs and supernova of type Ic (i.e. core-collapse but no 
hydrogen and helium lines). These type of stars are usu- 
ally identified as Wolf-Rayet and are characterised by initial 
masses of about 25 - 40 M Q (Maeder & Conti 1994). The 
heavy mass loss rate that distinguishes these objects reduces 
the mass of the star by losing the outer hydrogen and he- 
lium layers, leaving a carbon/oxygen core (Matheson 2004). 
In general they are expected to have less massive envelopes 
than SNe lb. Typical masses are around 16 to 18 Mq but 
the range is wide: from 5 Mq to 48 Mq . Therefore W-R 
stars are the low mass descendants of previously massive O 
stars (Maeder 1990). 

Radii are very difficult to determine for these stars, es- 
pecially because of the strong mass loss that makes it diffi- 
cult to define what is the stellar surface. However, few esti- 
mates exist for eclipsing binaries, where the radius is found 
to be 11 R Q for the late type CQ Cephei (HD 214419, 
WN7+09) and ~ 37? for the earlier V444 Cygni (HD 
193576, WN5+06). There appears to be a correlation be- 
tween large radii and late-type WR stars (where the Hy- 
drogen lines are detected), while early type WR stars (no 
Hydrogen lines) show smaller radii. 

Therefore such a smaller progenitor with a thinner stel- 
lar envelope would have more favorable conditions for the jet 
to emerge from the stellar surface (as for the GRB 021004, 
see Starling et al. 2005). GRB 031203, GRB 030329, and the 
most popular GRB event related to SN 1998bw, all show 
spectra typical of a type Ic SN. 

Another possibility would be that the jet is already 
present before the SN event, as one can see in young stellar 
sources, whose jet are well observable and often in preces- 
sion, such as Herbig Hero objects. 

Although there is no generally accepted mechanism able 
to explain the making of a thin jet structure, there seems to 
be a link between jets and accretion discs. 



3.2 Muon jets to overcome the SN-GRB opacity 

In our previous models of GRB emission we assumed a jet 
made of electron pairs of energy around 1-10 GeV produc- 
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Figure 4. The electron and muon interaction lengths. The 
dashed- dotted and dotted lines correspond to the synchrotron en- 
ergy loss distance (for muons and electrons respectively) for dif- 
ferent values of the magnetic field: 100 G, 1 G and 10 — 2 G. The 
straight solid line labelled t^ indicates the muon lifetime. The 
dashed lines indicate the IC interaction lengths for muons and 
electrons. Finally the two solid curves labelled /i + fi~ and e+e — 
correspond to the attenuation length of high energy photons pro- 
ducing lepton pairs (either ^ or e^°) through the interaction 
with the SN radiation field. We have assumed that the thermal 
photons emitted by the star in a pre-SN phase have a black body 
distribution with a temperature T ~ 10 5 K. Assuming a radius 
R ~ lO-f?0, we are considering a luminosity of Lsn — 2.5 ■ 10 41 
erg s — 1 . Around 10 15 — 10 16 eV muons decay before losing energy 
via IC scattering with the stellar background or via synchrotron 
radiation. 



10 -3 s (~ 300 km) from the stellar surface the baryon den- 
sity would also be able to reduce and block the UHE electron 
pair propagation. This represents a strong constraint for a 
jet model made of e . 

UHE muons (E^ > PeV) instead are characterised by a 
longer interaction length either with the circum-stellar mat- 
ter and the radiation field, since the fi — 7 cross section 
scales as the inverse square of the mass of the lepton involved 
(in the Thomson regime). In the Klein-Nishina regime at E 
> 10 17 eV the interaction lengths for muons and electron 
becomes comparable (see Fig. ■ 

For this reason here we discuss the possibility that the 
primary particles of the jet may be muons with _E M ~ 10 lj 
- 10 16 eV, which can more easily escape from the stellar in- 
terior. A jet of primary electrons at the same energy would 
have a similar IC interaction length but their escape from the 
stellar surface would be obstructed by both the baryon load 
(i.e. electromagnetic interactions inside the matter) and by 
the synchrotron opacity which would be dominant at those 
energies for a magnetic field between 1 and 100 Gauss (see 
dotted lines in Fig. From this starting point we investi- 
gate the mechanism able to produce the observed gamma 
radiation in GRBs from muon jets. 

We assume that the star where the GRB occurs is a 
WR-like object with a characteristic radius R t ~ WRq. 
When the muon jet breaks through the surface of the star 
it may interact either with the magnetic field of the star or 
inverse Compton (IC) scatter its optical - UV photon field. 

The importance of the synchrotron process depends on 
the magnetic field configuration around the jet and the star. 
In the case of a WR star, the magnetic field on the surface 
may be as high as 100 G and assuming a 4j dependence, 
one can estimate the field intensity at a distance r. 

As an example, for an initial muon at 2? M ~ 10 16 eV and 
a magnetic field of roughly 1 Gauss, synchrotron photons 
would be emitted at 



ing the gamma-ray emission by IC onto thermal photons 
(Fargion 1995a, 1995b, 1995c; Fargion 1998; Fargion 1999; 
Fargion 2003). However, the main difficulty for a jet of GeV 
electrons is that their propagation through the SN radiation 
field is highly suppressed. 

Indeed high energy electrons and positrons (E e > 1 
GeV) may inverse Compton scatter the UV/optical photons 
emitted by the star to high energies leading to an electron- 
photon cascade. The interaction length of e ± of energy E e , 
in an isotropic photon gas with density n(e), is given by 
(Protheroe 1986) 



A 



ic 



n(e) 



c(w)(l — /3 cos 6)d(cos 9)de 



(2) 



where 9 is the angle between the electron and photon 
directions in the laboratory frame and lo = (tE /m e c 2 )(l — 
(5 cos 9) is the photon energy in the electron rest frame. 

In Fig. [I] we show the IC interaction length for elec- 
trons and muons [dashed lines) scattering a black-body radi- 
ation field at a temperature T = 10 5 K. In the Klein-Nishina 
limit we have used the asymptotic formula given by Gould 
& Raphaeli (1978). From Fig.|I]one can see that electrons at 
E e > 10 GeV can not propagate for more than 10 -3 s due 
to the scattering with the optical/UV emission. Moreover at 



rpsync 
-C/ 7 



4 x 10 



5 f- 



E u 



(re) 



eV 



VlO 16 eV , 

and the characteristic interaction length is 



E 



(4) 



It follows (see Fig. that for B < 100 G the syn- 
chrotron interaction distance scale (dashed- dotted lines) is 
much longer than the IC (upper dashed line) and the muon 
decay length (straight solid line), therefore muons can prop- 
agate outside the star without significant synchrotron losses. 
Only for a magnetic field around 10 3 G, muons would lose 
most of their energy to synchrotron emission of photons with 
E 1 ~ 400 MeV, but such high magnetic fields in stars are 
yet to be confirmed. 

Nevertheless for E^ ~ 10 16 eV even the interaction 
length of the IC scattering onto the stellar radiation field is 
larger than the muon lifetime (Fig. , therefore such muons 
decay into electrons (in roughly 100 s) before their energy is 
dissipated via IC or synchrotron. As a consequence, a jet of 
muons would be able to escape from the star and propagate 
for about 100 light-seconds before decaying into electrons. 
To understand the fate of such high e ± , one needs to com- 
pare again the IC and synchrotron interactions to establish 
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which is the mechanism that most affects their loss of energy 
as they propagate outside the stellar surface. 



4 GRBS BY ELECTRONS SHOWERING IN 7 
JETS 

The main source of the thin gamma jet is the ultra high en- 
ergy electron pairs showering into thin 7 jets. PeV electrons 
and positrons from muon decay may either inverse Compton 
scatter the UV/optical photons emitted by the star, leading 
to an electron-photon cascade, or lose energy emitting syn- 
chrotron radiation in galactic magnetic fields. 

From Fig.^Jit appears that synchrotron interactions are 
less favoured compared to ICS because their distance scale is 
larger. In fact, even assuming that intensity of the magnetic 
field of the pre-SN star, B, decreases as r -2 (rather than 
r -3 ), at 100 light seconds from the star, B is low enough 
(< 10 -2 G, if the field on the surface of the star is ~ 100 
G) that electrons survive synchrotron losses, and IC consti- 
tutes the dominant energy loss mechanism (see in Fig. 0]the 
synchrotron energy loss distance for B — 10~ 2 G and the IC 
interaction length for electrons at energies around 10 eV). 

Compton interactions for PeV electrons occur in the 
Klein-Nishina regime, since the ambient photon energy in 
the rest frame of the electron is u > m e c 2 , and the e 
pair tend to transfer most of their energy (about \E e ) to 
the background photons. The propagation of the hard 7- 
rays created is attenuated by the collisions with the stellar 
radiation background that leads to pair production (7 + 7 — > 
e + e~). Compton scattering and pair production interactions 
concur to produce a complicated "electromagnetic cascade" 
which appears to be the favoured mechanism for the PeV 
electrons to lose most of their energy (see also Fargion & 
Colaiuda 2004). 

The attenuation length for a hard photon of energy E, 
traversing a black-body photon gas of energy e and density 
n(e) is given by a relation similar to Eq. [5] 

X pair = ^J n(e) J cr 77 (s)(l -/3cos6>)d(cosl9)de (5) 

where cr(s) is the pair production cross section for pho- 
tons with center of mass energy %/s, and s = 2eE(l — f3 cos 6). 
Apair has been derived following Brown et al. (1973) and it 
is shown in Fig. [I] (solid curve labelled e + e~) as well as the 
interaction length of the process 77 — > . 

To summarise, in the e — 7 collisions (ICS), PeV elec- 
trons transfer half of their energy to the ambient pho- 
ton, producing a hard 7-ray which then creates an e + e~ 
pair due to the interactions with the photon stellar back- 
ground (7 — 7 — > e + e~). However, as the electron propa- 
gates outward, the photons encountered are progressively 
redshifted, and the interaction length of the e 7 and 77 
scattering is continuously changing with distance and time 
as we show in Fig. [K] We are assuming that the energy of 
the photons is redshifted by a factor (1 — /3cos9), where 
tanf? = Ri,/(l e + R*), with l e being the distance covered 
by the electron at a certain time t, and 7?* is the WR star 
radius. The multiple photon-electron cascade continues dur- 
ing the electron propagation. For an initial PeV electron, we 
consider that about 10 collisions are needed for its energy 




Figure 5. The evolution of the e-y (dashed line), 77 (solid line) 
interaction length as a function of time. Here we have plotted 
the curve at t = s, 100 s (when ^ decay into e^), 200 s (the 
lapse time of ten e7 bounces), and 1000 s (from left to right). The 
IC and pair production interaction lengths are compared to the 
synchrotron one for B = 100, 1, 10~ 2 , 10" 4 Gauss. 



to be degraded to E e ~ 500 GeV - 1 TeV (see also Far- 
gion & Colaiuda 2004). The lapse time for ten bounces is 
approximately 10 x 10 s ~ 100 s. At this energy, the photon 
produced at E 1 = \E e is able to "escape" without inter- 
acting anymore with the stellar background (Fig. |SJ . This 
mechanism would imply a hard 7 emission (TeV, sub- TeV) 
from GRBs. 

Seven GRBs with emission beyond 100 MeV have been 
detected by EGRET on board of the Compton Gamma Ray 
Observer (CGRO) (Schneid et al. 1992), including the long- 
duration burst (~ 5000 s) GRB 940217 whose peak emis- 
sion has been measured at ~ 18 GeV (Hurley et al. 1994). 
In principle one cannot exclude that GRBs may emit at 
TeV energies, even though the propagation of these photons 
would be suppressed by the presence of the infrared (IR) 
extragalactic background, (jTeV + Jir — * e + e~), so that 
only photons from relatively nearby sources (< 100 — 300 
Mpc) could be detected. A few claims of GRB detections in 
the TeV and sub TeV range have appeared in the literature 
in the last few years. Such high energy signals have been 
claimed to be observed by Milagrito in one out of 54 GRBs 
from the BATSE catalogue (Atkins et al. 2000). Another 
claim of possible sub- TeV emission from GRB 971110 has 
come from the GRAND project (Fragile et al. 2004). How- 
ever one has to keep in mind that all these results are not 
considered as firm detections. It follows that even if we may 
expect a very high energy (~ TeV) "precursor", it would be 
very difficult to observe it because of the IR cut-off. 

Going back to the keV - MeV emission, in this scenario 
it would result from the continuous IC scattering of the same 
electrons, when they have reached lower energies. In fact, af- 
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ter they have been degraded to few hundreds GeV, the elec- 
trons keep interacting with the stellar background but the 
process is not anymore in the Klein-Nishina regime. These 
collisions lead to the creation of photons of energy given by 



fif = ^ 7 2 e(l-/3cos(9) 



(6) 



where is the scattering angle between the photon and 
the electron 3 . For tg6 = R+/(l e +R*) < 1, (1-/3 cos 6) ~ 6 2 
and the equation reduces to 

where the 6 2 factor accounts for the apparent redshift 
in the electron frame. The energy loss rate in this case has 
the form 



dt \ 3 m|c 4 



-bE z 



(8) 



and the energy decrease as a function of time is given 



by 

E(t) = {E(t y 1 + bty\ 



(9) 



Thus the collisions with the ambient photons rapidly 
lower the energy of the electrons. When their energy is re- 
duced to E e = WGeV, from Eq.0 one obtains approximate 
photon energies of 450 KeV (j^y) 2 (^±y) ( 10 Xad Y- 



5 THE SGR 1806-20 GIANT FLARE BY A 
BLAZING GAMMA JET 

An isotropic release of energy, as predicted by the magne- 
tar model, shows some inconsistencies when applied to the 
properties of the giant flare from SGR 1806-20, as we have 
pointed out in §1.3. There are two more issues related to 
the magnetar interpretation that we want to discuss in this 
section. 

First, how can a spherically symmetric explosion jus- 
tify the peculiar features observed in its light curve? A pre- 
cursor burst occurred 143 s before the initial onset, lasting 
slightly longer than 1 s, and it was followed by a very in- 
tense spike with a very short timescale (0.7 sec), the blazing 
giant flare that saturated most GRB detectors. After t = 
61 until 170 ms the signal " decreased gradually with oscil- 
latory modulation, which suggests repeated energy injections 
at ~ 60 ras intervals" (Terasawa et al. 2005). Soon after 
the peak of the emission, between 400 and 500 ms, several 
humps were detected in the decay profile from both Swift 
(Palmer et al. 2005) and GEOTAIL (Terasawa et al. 2005). 
Then the profile declines again until t ~ 300 s, with the light 
curve showing signs of pulsations with a 7.56 s period. An 
increase in the emission seem to occur again around t ~ 400 
s, with a long bump peaking at t ~ 600 - 800 s (Mereghetti 
et al. 2005). An even longer tail seems to appear from the 
data, which lasted until 1000 - 2000 sec after the first trigger 
(Mereghetti et al. 2005). 

Secondly, does the relation between PP and B 2 , pre- 
dicted by the magnetar model, find a confirmation from the 

3 A similar showering process has been discussed for UHECR 
from BL Lac sources (Fargion & Colaiuda 2004) 



observations? The magnetar scenario assumes a neutron star 
with a huge magnetic field that is estimated by observing the 
rotational energy loss (mainly due to the dipole moment ra- 
diation) and the rate of rotation variability. 

Let us calculate these quantities for SGR 1806-20 as- 
suming a simple dipole model. The rotational energy and 
its derivative (for a nominal period of Psgr = 7.56s) are: 



E rot = ^Inslo 2 ~ 3.6 • \Q 4i p- 2 



Ins 



10 45 <? cm? 



E rot = Insuu = 4 • 10 46 PP" 3 



V 10 45 # cm 2 



erg 



erg s 



(10) 



(11) 



However from the dipole radiation formula one finds: 



Erot — ~ 



{2vfB 2 NS R%, 
6c 3 P 4 



sin 2 (0) 



(12) 



From the last two equations one derives that the mag- 
netic field is related to the period P and its derivative P as 
follows: 



B = 3.2 • 10 19 (PP)5 G 



(13) 



Mereghetti et al. (2005) have shown the long term P 
variation of SGR 1806-20 between 1993 and October 2004, 
using data from ASCA, RXTE, Beppo Sax and XMM. The 
period derivative increased from 8x 10 -11 s s~ x (before 1998; 



Kouveliotou et al. 1998) to (5.49 ± 0.09) x 10~ lu s s 



- 1 (be- 



tween 1999 and 2004), while the period P changed from 7.47 
s, to 7.56 s. This implies that the magnetic field has changed 
from 



Bin 
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to 
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G (15) 



it follows that PP has increased roughly of a factor 7, 
as the square of the magnetic field, B 2 . According to these 
values, the corresponding total energy E ~ ^B 2 R% S must 
have changed from 



E m = 1.0-RT 



to 



E fin = 7.1-10 4 



(7.47 s) 
(7.56 s) 



8- 10- 



(wt) -»w> 
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Where does such an increase (7 times higher) in the 
energy density come from, given that the rotational energy 
is at least three orders of magnitude lower (see Eg. HOB ? 

Let us go back to the intense flare activity of 2004 
December 27. No sudden change in the spin frequency 
(y = 1/P) was found after the flare and an upper limit 
of Au < 2 x 10~ J Hz has been set on the frequency varia- 
tion (Woods et al. 2005). On the other hand a decrease of a 
factor 2.7 in the spin-down rate has been recently observed 
roughly one month after the giant flare, compared to the 
average value of the past 4 years (Woods et al. 2005). If PP 
is proportional to B 2 , the decrease in P should imply from 
Eq. H7l that the energy released in the flare should have been 
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Figure 6. A close up of the 2D trajectory of the processing jet 
and of the way it blazes the observer along the line of sight (upper 
panel ), and the corresponding light curve profile (lower panel ). 
Here we have used the same set of parameters described in §2 
and in Fig. [5] This set provides a good fit to the light curve 
of SGR 1806-20 observed during the giant flare (Hurley et al. 
2005, Mereghetti et al. 2005). The early precursor detected ~ 100 
s before the onset of the giant flare, may be due to a random 
trajectory of the jet not shown in this picture. 



~ 5 x 10 47 erg, while only a tenth of this value (~ 5 x 10 46 
erg) has been recorded by RHESSI (Hurley et al. 2005) and 
GEOTAIL (Terasawa et al. 2005) . Where has the remaining 
energy gone? 

It follows that the evidence for a correlation between 
PP and B 2 is not at all clear thus far, especially within the 
magnetar scenario. 

Again, we suggest that a thin collimated jet is able to 
solve these puzzles better than the magnetar model. 

In Fig.|Sland|7| we show the results of our model (see Eq. 
with the set of parameters displayed in Tabled First, a 
thin jet would be able to lower the total energy budget of the 
process {Esgr ^ sn-^W ~ 10 9 x 10 38 erg s" 1 ). Secondly, 
a jet in precession would provide a natural interpretation of 
the features observed in the light curve of SGR 1806-20 as 
one can see from Fig. |S| 

The structure of the 7 jet (originated either via ICS 
or synchrotron radiation) consists, because of the relativis- 
tic kinematic, in concentric cones where the inner and more 
collimated ones correspond to the more energetic gamma 



Figure 7. A simulation of the SGR 1806-20 activity along a time 
range of ~ 30 yr showing a periodicity of ~ 6.5 yr. This result has 
been obtained with the same set of parameters discussed in §5.1. 
The clustering of the events is related to the precessing motion of 
the entire cone described by the jet, whose period is 2ir/ui a 13 

yr- 



radiation. A section of this cone perpendicular to the axis 
would show different concentric rings. If the motion of the jet 
is slow 4 , the profile of the light curve would appear symmet- 
ric to an observer, whereas if the jet is fast, the inner struc- 
tures (at higher energy) move quicker than the outer ones 
(at lower energy). This determines a "compression" of the 
high energy cones along the direction of the precessing mo- 
tion, the faster the motion of the jet the stronger this effect. 
Fast moving jets though should appear with a fast-rising 
hard signal, followed by a lower energy tail. Consequently, 
a fast thin jet that suddenly blazes the observer along its 
precessing motion may explain the short e-fold time of the 
initial rise (At ~ 1.3 ms; Terasawa et al. 2005). The sudden 
oscillatory modulation detected along the decay profile at 
61 ms < t < 170 ms and 397 ms < t < 500 ms (Terasawa et 
al. 2005; Palmer et al. 2005) finds explanation if we assume 
At ~ P/7 ~ 1.3 • 10 -3 s for P = 7.5 s and 7 ~ 10 4 . 

We imagine that the rotational energy losses of the pul- 
sar are mainly related to the braking torques of an accretion 
disc which determine the change of the period derivative. 
This implies that PP is not correlated to the magnetic en- 
ergy density (contrary to the magnetar model) before and 
after the flare. 

The recent increase of the hardening might be due to 
the motion of the cone described by the precessing jet: the 
more it bends towards the observer the longer is its P, the 
more frequent may be the SGR activity and the harder is 
the spectum (Mereghetti et al. 2005). 

Moreover the observations of the radio afterglow of the 
SGR 1806-20 flare (Taylor et al. 2005) give evidence for a 
polarization variability, for a rebrightening in the radio light 
curve (Gelfand et al. 2005) and they show a variation in the 
position of the centroid of the afterglow, that appears as 
asymmetric and elongated: all these properties are in favour 



4 Fast and slow jets are defined by the ratio v±/c (> or < 1 
respectively) where v±_ = u] psr RL svci6 psr and Rl ~ c/P 
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of the jet interpretation. In this case a wider and less colli- 
mated cone of GeV elecrons may produce the radio emission 
of the afterglow. 

In our model, highly energetic particles are ejected from 
the poles of a spinning neutron star. Here we propose three 
possible scenarios that can account for the gamma emission 
from SGRs, particularly for what has been observed in the 
7-ray giant flare from SGR 1806-20. (1) In the simplest case 
the jet is made of GeV electrons that emit 7 rays via ICS 
onto thermal photons, as we have already discussed in pre- 
vious works. Then we examine the possibility that the soft 
gamma flare results from synchrotron radiation (due to the 
galactic magnetic field) of PeV electrons possibly originated 
by either (2) a primary jet made of muons, or (3) by EeV 
nucleons (p + 7 — > TT + n;n ^ /j, ^ e; n ^ p + e + i? e ). At this 
stage we cannot decide which model is able to best interpret 
the giant flare emission. PeV muons have the advantage of 
escaping for about 100 seconds before they decay. On the 
other hand the introduction of an hadronic component of 
the jet might be necessary, given the evidence for a spatial 
correlation between an excess of EeV CRs and the location 
of SGR 1806-20. 



5.1 Gamma jets by IC of GeV electrons 

As we proposed in previous works, the simplest approach 
is to assume that the gamma emission arises from Inverse 
Compton Scattering of GeVs electron pairs (j e 2 x 10 3 ) 
onto thermal photons (Fargion 1995, 1996, 1998, 1999). 
Their ICS will induce an inner jet whose angle is A8 < I/7 ~ 
5- 10~ 4 rad ~ 0.0285° and a wider, less collimated X, optical 
cone. Indeed the electron pair Jet may generate a secondary 
beamed synchrotron radiation component at radio energies, 
in analogy to the behaviour of BL Lac blazars whose hard- 
est TeV component is made by ICS, while its correlated X 
emission is due to the synchrotron component. Anyway the 
inner jet is dominated by harder photons while the external 
cone contains softer X, optical and radio waves. In a first 
approximation the gamma emission is given by the Inverse 
Compton relation: < e 7 >~ -y^kT for kT ~ 10" 3 - 10" 1 
eV, and E e ~ GeV, leading to the characteristic X spec- 
trum. However the main puzzle of this scenario is the way 
of originating such high energy particles and how they can 
escape from the strong magnetic field of the pulsar. Other 
mechanisms and progenitors need to be introduced to guar- 
antee the escape from the star of the high energy electrons 
forming the jet. 



5.2 PeV Muon jets in SGRs ? 

Similarly to the GRB case, we assume a primary jet made of 
muons with £7 M ~ 1 - 10 PeV, that can travel for about 100 
s before they decay into electrons (see Fig. ■ The presence 
of a UV/optical photon background that can prevent the 
propagation of high energy electrons is negligible around a 
pulsar. 

We are also assuming that the jet is aligned with the 
magnetic lines, thus muons would not emit synchrotron pho- 
tons as they move through the very intense magnetic field 
of the pulsar generally expected near its surface (possibly as 
high as 10 15 G). The advantage of a primary beam of muons 



is that it guarantees that the electrons are produced away 
from the surface of the neutron star, where the magnetic 
field would be much lower. 

When electrons of energy ~ IPeV arise from the de- 
cay of muons after about 100 s, the pulsar magnetic field 
at this distance has decreased by a factor Bi^{R psr /j*)~ 3 
G ~ 3 x 10- 5 B 15 (R par /10km)- 3 G, with B15 = (B/10 15 
G). In practice B at 100 s is comparable to the galac- 
tic magnetic field. Given that the synchrotron losses for 
B ~ 10~ 5 — 10~ 6 G acts on a fairly long time scale, it is 
more likely that PeV electrons lose most of their energy 
while they propagate through the galaxy and they interact 
with the galactic magnetic field. This would lead to a radi- 
ation of energy 

with an energy loss distance scale of 

The Larmor radius for the galactic magnetic field is 

and it is about three orders of magnitude lower than the 
synchrotron energy loss distance (RL/ct sync = 1.3 • 10 -3 ). 
This may cause a widening of the opening angle of the jet, 
that would have a final "fan"-like (see also Fargion 2002, Far- 
gion et al. 2004) or a double disc where the e are bounded 
by the magnetic field. The jet would open along a two di- 
mensional plane, given that the galactic magnetic field is 
homogeneous on a scale of 55 pc (Protheroe 1990). We are 
assuming a Lorentz factor which is j e ~ 2 x 10 9 , but in this 
case the solid angle of the jet is not oc I/7 2 ~ 2.5 x 10 -19 , 
but I/7 x 9l, where 9l is the aperture angle due to the 
Larmor precession around the galactic field. For Ql — 10°, 
An ~ 10 °, and this aperture angle guarantees that the 
number of events is not suppressed by a too narrow beam- 
ing. 



5.3 EeV nucleons as progenitors of SGR 7 jet? 

Finally we consider the possibility that SGRs may be sources 
of UHECRs, and that the PeV muons may be originated 
by jets made of higher energy hadrons. Similar ideas on 
UHECRs have been recently suggested also by other au- 
thors (Asano, Yamazaki, & Sugiyama 2005), but here we 
assume that the energy of the hadrons from SGRs does not 
exceed 10 18 eV instead of being around 10 20 eV as proposed 
by Asano and collaborators. Indeed there is strong evidence 
in CR spectra that the origin of EeV particles may be both 
galactic and extragalactic. Excluding very rare SN/GRBs, 
in our galaxy there are no active sources at this range of 
energies, thus we cannot ignore the role of SGRs as possibly 
related to EeV CRs. 

In 1999 AGASA reported the measure of a significant 
anisotropy (4%) in the arrival direction of CRs in the small 
energy range 10 17 ' 9 — 10 18 ' 3 eV as due to the excess of 
events in two regions of about 20° near the SGR 1806-20 
and Cygnus, and it has been argued that this anisotropy 
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Figure 8. The correlation between the BATSE data, the AGASA 
discovery of an anisotropy in the arrival direction of EeV CRs 
near the galactic center and the Cygnus regions and the position 
of the galactic SGRs. The four clusters of the data from BATSE in 
the left-hand side of the map correspond from top to bottom to: 
the Cygnus region, SGR 1900+14, SGR 1806-20, SGR 1627-41. 
Note that AGASA was unable to record of the last SGR because 
below its horizons. The fifth SGR, 0526-66 in the Magellanic cloud 
is not indicated in this map because it did not show signs of 
activity between 1997 and 2000. Most of the authors Hayashida 
et al. (1999) or Fargion (2002) argued that this excess in the CR 
distribution might be due to EeV neutrons partially decaying in 
flight. 



was due to EeV neutrons (Hayashida et al. 1999). Interest- 
ingly enough we have also found a correlation between the 
BATSE data of SGR 1806-20 and the AGASA excess, as one 
can see from Fig. |H1 where we have overlaid the BATSE de- 
tections (between 1997 and 2000) to the AGASA map. This 
would suggest that the neutron signal and the soft-gamma 
detection from BATSE are related to the same mechanism. 

To explain such an interesting correlation, we will start 
with the assumption that protons can be accelerated by 
the large magnetic field of the pulsar up to EeV energy. 
The protons could emit directly soft gamma rays via syn- 
chrotron radiation with the galactic magnetic field (E? ~ 
10(E p /EeV) 2 (B/2.5 ■ 10~ 6 G) keV), but the efficiency is 
poor because of the too long timescale of proton synchrotron 
interactions. Photopion production must occur to produce 
the observed neutron excess, a process that also gives birth 
to neutral and charged pions. Which photons represent the 
target for such a process is not completely clear yet. It is 
possible that the object may have an IR emission and that 
this thermal background acts as a target for the high en- 
ergy protons. In the literature there are claims of identifi- 
cation of an IR counterpart to SGR 1806-20. Eikenberry et 
al. (2001) report the discovery of possible IR counterparts 
to the object, and more recently Israel et al. (2005) claim 
the detection of one faint object (K 3 = 21.6 ± TO) which 
is compatible with the radio position of the SGR (Gaensler 
et al. 2005) obtained with a very accurate positional accu- 
racy (of 0.1"). The energy of the thermal photons neces- 
sary to produce pions is e 7 = 0.2 GeV 2 /E p = 0.2 eV for 
E p = 10 18 eV. Charged pions (born with roughly a tenth 
of the original energy of the proton) decay into muons and 



then into electrons with E e < 10 eV. Medina- Tanco & 
Watson (2001), following Stecker (1968) argues that dur- 
ing the propagation of EeV protons from the galactic cen- 
tre, photopion production with ambient galactic IR photons 
(p + Jir — > A + — ► n + 7r + ) is favoured compared to proton 
proton collisions (p + p^n + p + Ntv) in order to produce 
charged pions, and they obtain a mean free path of 37 pc 
for the interaction with an IR galactic photon background 
(Tir ~ 100 K). Thus, even if there is not an IR source lo- 
calised with the SGR, the interaction with the galactic IR 
background may be efficient in converting the protons into 
pions and neutrons (see also Grasso & Maccione 2005). 

Pions then decay into muons that decay into high en- 
ergy electrons (E e <J 10 PeV) that, interacting with the 
local galactic magnetic field, lose energy via synchrotron ra- 
diation: 

- 4.2 x 10- (^ W ) 2 (m^—) eV (21) 
with an interaction length given by 

tavnc - ^ x 1010 ( 5-ifw) - 1 (2tw*g) 2 s (22) 

The Larmor radius is about two orders of magnitude 
smaller than the synchrotron interaction length and this may 
imply that the aperture of the jet is spread by the magnetic 
field. 

In this "hadronic" scenario we are assuming a Lorentz 
factor which is 7 e ~ x 10 10 , and, as mentioned in the pre- 
vious section, the solid angle of the jet is not oc I/7 2 ~ 
xlO -20 , but I/7 x 8l, where 0l is the aperture angle due 
to the Larmor precession around the galactic field. For 
6l — 10°, AS! ~ 10 -11 , and this aperture angle again guar- 
antees that the number of events is not suppressed by a too 
narrow beaming. 

The photopion interactions produce also neutrons at 
~ 10 18 eV (7 ~ 10 9 ), which have a decay length of < 
10 kpc (E/10 18 eV), and they can easily propagate along 
these distances without being bent by the galactic magnetic 
field. Neutron beta decay produces a proton, electron and 
an electronic antineutrino. Given its long decay length, the 
production of the electron would occur well far away from 
the original pulsar. The electron is produced with an energy 
spectrum that ranges between 10 15 eV < E e < 10 17 ' 5 eV and 
peaks at E e ~ 10 15 eV. A gain PeV electrons interact with 
the local galactic magnetic field producing via synchrotron 
radiation photons of few hundreds keV as given by Eq. 1181 

For higher electron energies, around E e > 10 16 5 eV the 
Larmor radius becomes comparable and even larger than the 
synchrotron interaction length, therefore the opening angle 
of the jet would remain collimated. However such an angle 
would be too small and would make SGR events extremely 
unlikely. Again for the higher energy electrons the peak en- 
ergy of synchrotron photons would be 

E r^l-^^(^)\^ 7 ^)eV (24) 

and to explain the hundreds keV emission we have to as- 
sume that the jet is not pointing exactly towards the Earth, 
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but it is slightly off-axis so that only the softer range of 
frequencies could have been observed during the flare. 



6 CONCLUSIONS 

The most popular models to date to interpret GRBs and 
SGRs, the fireball collimated in a few degrees jet for GRBs, 
and the spherically symmetric magnetars for SGR, can- 
not completely explain all the issues that these events do 
present. 

We have proposed an alternative scenario by introduc- 
ing a model with highly collimated jets of high energy parti- 
cles in precession. We are assuming the energy equipartition 
between the power of the jet and that of the SN (for GRB) 
or the X-ray pulsar (for SGRs). The gamma rays emerge 
in a small opening angle, and the narrow beaming allows to 
dramatically reduce the total output power (from ~ 10 50 erg 
s" 1 for the fireball/jet model of GRBs to 10 43 - fO 44 erg s _1 
in our GRB case as the typical SN power, and from 10 47 erg 
s _1 for the magnetar to fO 38 erg s _1 in our SGR scenario, 
comparable the the observed power of X-ray pulsars) , which 
is the most radical assumption in such isotropic models. 

Compared to our previous versions here we have con- 
sidered that high energy electrons are the decay product of 
a primary particle that in the case of the GRB is a muon. 
For SGRs we have considered a jet of either prompt muons 
or of protons that produce, through photopion production, 
secondary muons and electrons. This proposal arises from 
the fact that the propagation of a jet of prompt electron 
pairs would be severely suppressed by the extreme opacity 
conditions in the proximity of the surface of a pre-SN star 
or a highly magnetised pulsar. 

If GRBs are originated in a core collapse of very massive 
stars, a jet of muons would have a higher chance to escape 
from the stellar surface. Secondary electrons produced af- 
ter the decay of the muons at about fOO s from the surface, 
may produce the gamma emission after a long chain of inter- 
actions with the stellar background radiation (ICS), which 
gradually reduce their initial energy, from ~ 10 lj eV to a 
tens of GeV. A narrow jet would blaze the observer at high 
energy when viewed on-axis (at viewing angle ~ I/7), while 
it would appear "softer" if observed off-axis. This could give 
an explanation to the wide range of properties and emission 
from different GRB events. A larger viewing angle would be 
associated to the softer emission of the X-ray afterglow, and 
it may also be linked to the properties of objects such as the 
XRFs. 

A similar interpretation has been introduced for SGRs, 
with a particular attention to the giant flare from SGR 1806- 
20 on 2004 December 27. We have discussed different scenar- 
ios at the origin of this event. We have shown how a primary 
jet of muons (with _E M ~ 10 15 — 10 16 eV) decaying into high 
energy electrons, may be a source of a collimated gamma ra- 
diation. In this case the favoured mechanism of the emission 
would be the synchrotron interactions with the galactic mag- 
netic field, since the secondary electrons would be produced 
at a distance from the pulsar where its magnetic field would 
be negligible. The interactions with the galactic magnetic 
field will widen the shape of the jet like a "fan" (see also 
Fargion 2002, Fargion et al. 2004), and its solid angle would 
still be very narrow AQ oc 1/7 ~ 10" 10 - 10 -11 . In this 



scenario the PeV leptons (jj, decaying into electron pairs) 
might also be the secondaries of EeV protons that generate 
charged pions and neutrons via photopion production with 
the galactic IR background (Medina- Tanco & Watson 2001; 
Grasso & Maccione 2005). We have considered that SGRs 
may be sources of EeV protons and neutrons, on the basis 
of the AGASA discovery of an anisotropy in the UHECR 
signal in this energy range (Hayashida et al. 1999). This 
excess seems to be related to the area where SGR 1806-20 
is located (Bellido et al. 2001) and we also found a cor- 
relation between, the AGASA excess, the position of SGR 
1806-20, SGR 1900+14, Cygnus and the BATSE public data 
that recorded the activity of these objects between 1997 and 
2000 (see Fig.®. 

We foresee that if the role of nucleons as primaries of 
the soft gamma emission of SGR 1806-20 is correct, the gi- 
ant flare might also be source of a prompt (and possibly 
repeating) 7-induced shower (made by 7 photons from the 
decay of PeV neutral pions) that may have been detected by 
Milagro in correspondence with the 2004 December 27 flare. 
We also expect that a rich component of the EeV neutrons 
(or protons from their decay) might appear in the AUGER 
or HIRES detectors, in rough coincidence with this event. 
Because of the delayed arrival time of protons, one should ex- 
pect also a long and persistent UHECR afterglow. Finally a 
signal of secondary muons at PeV energies, induced by high 
energy neutrinos from the SGR, might occur in Amanda. 

To conclude, we imagine that if the precessing jet model 
gives a correct interpretation of the properties of SGRs, SGR 
1806-20 will still be active in the next months and years. 
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